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ABSTRACT 

We combined the spectroscopic information from the 3D-HST survey with Herschel data to characterize the Ha dust 
attenuation properties of a sample of 79 main sequence star-forming galaxies at z ~ 1 in the GOODS-S field. The 
sample was selected In the far-IR, at A=100 and/or 160 /rm, and only includes galaxies with a secure Ha detection 
(S/N>3). From the low resolution 3D-HST spectra we measured the redshifts and the Ha fluxes for the whole 
sample (a factor of 1/1.2 was applied to the observed fluxes to remove the [Nil] contamination). The stellar masses 
(M*), infrared (L/ij) and UV luminosities {Luv) were derived from the SEDs by fitting multi-band data from GALEX 
near-UV to SPIRE 500 ^m. We estimated the continuum extinction Estor(B-V) from both the IRX=L/ij/L( 7 v ratio 
and the UV-slope, /3, and found an excellent agreement between the two. The nebular extinction was estimated from 
comparison of the observed SFR^fa and SFRjjv We obtained /=Esft,r(B-V)/E„e5(B-V)=0.93±0.06, i.e. higher than 
the canonical value of /=0.44 measured in the local Universe. Our derived dust correction produces good agreement 
between the Ha and IR+UV SFRs for galaxies with SFR> 20 M 0 /yr and M* > 5 x 10^° Mq, while objects with lower 
SFR and M* seem to require a smaller /-factor (i.e. higher Ha extinction correction). Our results then imply that the 
nebular extinction for our sample is comparable to that in the optical-UV continuum and suggest that the /-factor is 
a function of both M* and SFR, in agreement with previous studies. 

Key words, galaxies: high-redshift - galaxies: star-formation - dust, extinction 


1. Introduction 


The rate at which a galaxy converts gas into stars, the Star 
Formation Rate (SFR), is a fundamental quantity to char¬ 
acterize the evolutionary stage of the galaxy. The SFR can 
be measured in several ways, e.g. from the UV luminosity, 
the far-IR emission, the recombination lines or the ra dio- 
continuum ( Kennicutt|1998 Madau &: Dickinson]2014 1, al¬ 
though each of these SFR tracers suffers from some un¬ 
certainties. Concerning the UV and optical indicators, the 
main source of uncertainty in the estimate of the SFR 
is due to the dust extinction, that strongly absorbs the flux 
emitted by stars at UV and optical wavelengths and re¬ 


emits it in the far-IR. Hence, an accurate quantification of 
the impact of dust on the galaxy integrated emission is cru¬ 
cial to precisely evaluate the SFR. The dust distribution 
inside a gal axy can be described by a two compo¬ 
nent model dCharlot &: Fall 2000] e.g.), including a 
diffuse, optically thin component (the Inter-Stellar 
Medium, ISM) and an optically-thick one (the birth 
cloud) related to the star forming regions. The birth 
clouds have a finite lifetime {tbc ^ 10^ yr) and con¬ 
sist of an inner HII region ionized by young stars 
and bounded by an outer HI region. This model as¬ 
sumes that the stars are embedded in their birth 
clouds for some time and then disrupt them or mi- 
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grate away into the ambient ISM of the galaxy. In 
this model the emission lines are produced only in 
the HII regions of the birth clouds as the lifetime 
of the birth clouds is in general greater than the 
lifetimes of the stars producing most of the ioniz¬ 
ing photons (~ 3 X 10® yr). The emission lines and 
the non ionizing continuum from young stars are 
attenuated in the same way by dust in the outer HI 
envelopes of the birth clouds and the ambient ISM 
but, since the birth clouds have a finite lifetime, the 
non ionizing continuum radiation from stars that 
live longer than the b irth clouds is attenua ted only 
by the ambient ISM (Chariot &c Fall 2000). 


This two component model was conceived to explain the 
higher extinction observed on the nebular lines with respect 


(e-g- 

Fanelli et al. 

1988 

Mas-Hesse & Kunth[|1999 

Mayya 

et al.||2004[ jCid Fernandes et al.||2005[ [Calzetti et a 

l.[[1994[ 

2000[). 


Forster Schreiber et al. 2009 


[Whitaker et al.||2014[ |Yoshikawa et al.||2010[ tor samples 
of mainly optical and/or near-IR selected sources). Some 
recent works also suggested that the factor f=F,star{^- 
V)/E„eb(B-V) 

required to reconcile the various SFR measurements 
is higher than those computed in the local Universe (e.g. 
[Kashino et al.|p013| jPannella et al.|[2014[ on sBzk selec- 
tions). These contrasting results found in the literature 
are not surprising, given the different and indirect meth¬ 
ods used and/or the small and sometimes biased samples 
of most of these studies. Direct measurements of the 
Balmer Decrement would be required on statistical 
samples of sources at high redshift to clarify dust 
properties of distant star forming galaxies. 

In the present paper we used near-infrared spectros copic 
data from the 3D-HST survey (Brammer et al. 2012) and 
multi-wavelength photometry to derive the differential at¬ 
tenuation on a sample of 79 star forming galaxies at z be¬ 
tween 0.7 and 1.5, selected in the far-IR within the GOODS- 
South field. The presence of the far-IR photometry for the 
whole sample is very important, as these data allow us to 
robustly constrain the integrated dust emission. This kind 
of approach has been often appl ied to galaxies in the 
local Universe (see for example Dominguez Sanchez 


et al. [2014), but only very occasionally at higher 
redshifts. 

The paper is organized as follows. In Sectionthe prop¬ 
erties of the sample and the data-set are described. Section 
[^ presents the spectral analysis. Section [^ illustrates the 
computation of the physical quantities for the galaxies in 
the sample. Section describes the measurement of dust 
extinction on the Ha emission. Section [^ presents a critical 
discussion about the results whit carefully attention to the 
assumption behind our analysis and Section summarizes 
the results. 

We adopt throughout this work a standard cosmology 
(Ho = 70km/s/Mpc, = 0.3, Ha = 0.7) and assume a 
Salpeterj ( |I9^ IMF. 


The relation between the color excess of the nebular re¬ 
gions and that of the stellar continuum derived by [Calzettij 
erar] ( |T9M||20^ , i.e. f=E,tar{B-Y)/Eneb{B-V)=0.U, has 
proven to be successful for local star-forming galaxies, while 
it is still unclear if it holds true in the high redshift Uni¬ 
verse. 

The most secure method to quantify the amount of dust 
extinction in the HII regions is by directly measuring the 
Balmer decrement (i.e. the Ha/H/3 line ratio). 

Nevertheless, such measurements are very challenging 
at z > 0.5 since the Ha line is shifted to the less-accessible 
near-infrared window so indirect methods are often needed 
to infer the attenuation related to the nebular lines. 

Current studies on dust properties in high redshift 
galaxies lead to contrasting results. In fact, while some 
authors claim that the extinction related to the emission 
lines and to t he stellar continuum are comparab le in z ^ 2 
galaxies (e.g. Erb et al.|2006 Reddy et al.|20ld using UV- 
selection), other authors co iiiirm the validity of the Calzetti 
et al. local relation (e.g ~ 


2. Sample and data-set 

Our sample is selected in the far-IR using the Herschel 
obser vations from the PACS Evolutionary Probe survey 
(PEP, Lutz et al.|2011 1 in the GOODS-South field, that has 
the deepest sampling with PACS-Herschel photometry. The 
far-IR data are a key ingredient for the purpose of this work, 
as they allow us to strongly constrain the thermal emission 
by dust and therefore to infer robust estimates of the contin¬ 
uum dust attenuation for all the galaxies in the sample. The 
Herschel selected sample has been matched with objects at 
0.7 < z < 1.5 from the 3D-HS T survey catalog (Brammer 
et al.j2012[ [Skelton et al.|2014 ) for which the Ha line is ob¬ 
servable with the G141 grism. For the cross-correlation pro¬ 
cedure we took advantage of the GOODS-M ultiwavelength 
Southern Infrared Catalog {GOODS-MUSIC, Grazian et al. 
20061, that collects the available photometry from ~ 0.3 


to 8 fim for the objects detected in the GOODS-S field. In 
the following we briefly describe the catalogs used and the 
sample selection procedure. 

2.1. The Herschel/PEP survey 

The PEP survey is a deep extra-galactic survey based on 
the observations of the PACS instrument at 70, 100 and 160 
^m. The GOODS-South field is the deepest field analyzed 
by the PEP survey and it is the only one observed also at 
70 ^m. The 3(7 limit in the GOODS-S field is 1.0 mJy, 1.2 
mJy and 2.4 mJy at 70, 100 and 160 pm respectively. 

PACS photometry has been performed with a 
PSF fitting tool, adopting the positions of Spitzer 
MIPS 24 pm detected sources as priors. This ap¬ 
proach has been applied to maximize the depth of 
the extracted catalogs, to improve the deblending 
at longer wavelengths, and to optimize the band¬ 
merging of the Herschel photometry to the avail¬ 
able ancill ary data in the UV-to-nea r IR. As de¬ 
scribed in Berta et al. (2011| 2013[), the 24 pm 
have been used as a bridge to match Herschel to 
Spitzer/TR.AC (3.6 to 8.0 pm) and then to the opti¬ 
cal bands. PACS priors source extraction fol lowed 
the method described by Magnelli et al. (2009). The 
completeness of the catalogs have been estimated 
through extensive Monte-Carlo simulations and it 
turns to be on the order of 80%, with flux limits of 
1.39, 1.22 and 3.63 mJy at 70, 100 and 160 pm, re¬ 
spectively (see PEP full public data release]^. We 

http://www.mpe.mpg.de/ir/Research/PEP/DRl 
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defer to Lutz et al. (2011| for further information 
on the PACS source extraction performances. 


2.2. Multi-wavelength photometry 
2.2.1. The GOODS-MUSIC catalog 

The GOODS-MUSIC catalog reports photometric data for 
the sources detected in 2 : and Ks bands in the 
GOODS-South area and it is entirely based on pub¬ 
lic data. The catalog was produced with an accurate PSF 
matching for space and ground-based images of different 
resolutions and depths and it includes 14847 objects. 
The photometric catalog was cross-correlated with a master 
catalo^released by the ESO-GOODS’s team, that summa¬ 
rizes all the information about the spectroscopic red- 
shifts collected from several surveys in the GOODS 
area. For the sou rces lacking a spectroscopic measurement, 
(20061 computed a photometric red- 
technique with a set of 


Grazian et al 
shift using a standard 


synthetic templa tes drawn from the PEGASE2. 0 
synthesis model ( Fioc &: Rocca- Volmerange 11997). 

The MUSIG catalog was extended by Santini et al. . 


with the inclusion of the mid-infrared tiuxes, obtained from 
MIPS observations at 24 pm. This catalog was also used 
for the computatio n of the Spectral Energy Distributions 
(SED, see Section 4.1 for more details about the SED- 
fitting). 


2.2.2. GALEX, \RS-Spitzer and SPIRE observations 


In order to increase the spectral coverage of the final sam¬ 
ple, we added to the MUSIG photometry also the GALEX 
near-UV observations at A = 2310 A acquired from the on- 


the IRS observations at 16 pm ( 

Teplitz et al.|2005| and the 

SPIRE data at 

250, 350 and 500 pm (Oliver et al. |2012 

Roseboom et al.||2010HLevenson et al.||2010[ see also|Viero| 

et al.||2013 

and 

Gruppioni et al. 

20131. 


tremely important for the scope of this work. The GALEX 
data allow us to better constrain the SED at UV wave¬ 
lengths, enabling a more acc urate estimate of the UV spec¬ 
tral slope p (see Section 5.2 and Appendix [A] for details). 


On the other hand, the SPIRE data improve the spectral 
coverage in the far-IR, thus allowing us to derive a robust 
estimate of the bolometric infrared luminosity (see Section 


4 ^. 

The SPIRE catalog is also selected with posi¬ 
tional priors at 24 pm. SPIRE fluxes were obtained 
with t he technique described in [Roseboom et al.| 


(2010). Images at all Herschel wavelengths have un¬ 


dergone extractions using the same MIPS prior cat¬ 
alog positions. 

For the GALEX data we used publicly available 
photometric catalog from MAST and rely in this 
case on a simple positional association. 

The ancillary photometry was included by cross- 
correlating the above mentioned catalog using a 
positional association with a matching radius of 1 
arcsec. About 28 % of the galaxies in the final sample have 


^ www.eso.org/science/goods/spectroscopy/CDFS_Mastercat 
® http://archive.stsci.edu/index.html 


a near-UV observation, the 82 % of these sources have a 16 
pm counterpart, while 87 %, 62 % and 23 % of the galaxies 
have a SPIRE 250, 350 and 500 pm detection, respectively. 

This fast decreasing fraction with A is due to the de¬ 
grading of the PSF in the SPIRE diffraction-limited 
imager. 


2.3. The 3D-HST survey 


The 3D-HST n e ar-infrared spectros copic survey (Bram- 
mer et al. 2012 Skelton et al. 20141 covers roughly 75% 
of the area imaged by the GANDELS ultra-deep survey 


fields) Grogin et al. 2011, Koekemoer et al. 


I 


201 R. In this 


work we used the observations in the GOODS-S field from 
the preliminary data release, that covers an area of about 
100 square arcminutes. The WFG3/G141 grism is the pri¬ 
mary spectral element used in the survey. The G141 grism 
covers a wavelength range from 1.1 to 1.65 pm, thus can de¬ 
tect the Ha emission line in a redshift interval between 0.7 
and 1.5. The mean resolving power of this grism is about 
R ^ 130, insufficient to deblend the Ha and [NH] emissions. 
Our measured Ha flux is then affected by contam¬ 
ination from [NH], that we took into account and 
removed as described in Section 14.21 

Because no slit is used, and the length of the dispersed 
spectra is larger than the average separation of galax¬ 
ies down to the detection limit of the survey (Fobs (A) ~ 
2.3 X erg!sjcm? at 5cr), there is a significant chance 
that the spectra of nearby objects could be overlapped. 
This “contamination” by the neighbors must be carefully 
accounted for in the analysis of the grism spectra. To 
evaluate the contribution to the flux from nearby 
sources, we considered the quanti tative model de¬ 
veloped by Brammer et al. (|2012). For a complete 
description of the instrumental set-up, the data- 
reduction pipeline and the conta mination model 
we defer to Brammer et al. (2012|). Further details 
about the contamination issue are discussed in Sec¬ 
tions 12.41 and [s] 


2.4. Cross correlation and cleaning of the sample 

The use of the MUSIG catalog is fundamental for the asso¬ 
ciation between the Herschel and 3D-HST observations, as 
the Herschel’s beam is not directly matchable to the high 
resolution imaging of HST: the spatial resolution of PAGS 
at short and long wavelengths is ~ 5 arcsec and ^11 arcsec 
respectively, while the WFG3 camera has a spatial resolu¬ 
tion of 0.13 arcsec. The inclusion of the MIPS data at 24 
pm was very useful for the association with the Herschel 
data, due to the Spitzer - MIPS intermediate resolution 
between optical instruments and the Herschel Space Tele¬ 
scope. 

The PAGS-MUSIG catalog includes 591 sources, se¬ 
lected at 100 and/or 160 pm above the 3cr flux limits of 1.1 
mJy and 2 mjy, respectively. About 40% of these objects 
have ground-based spectroscopic redshift in the MUSIG 
catalog. We then cross-correlated the PAGS-MUSIG and 
3D-HST catalogs, using a matching radius of 2 arcsec, and 
found that 378 PAGS-MUSIG objects have a counterpart in 
the 3D-HST catalog. The choice of the adopted near¬ 
search radius, 2 arcsec, to cross-match the PACS- 
MUSIC sample and the 3D-HST is the best com- 
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promise between the 24 pm beam size (5.6 arcsec) 
and the IRAC ones (from 1.7 to 2 arcsec going from 
the 3.6 pm to the 8.0 pm channel). We have how¬ 
ever verified that using a slightly smaller (1 arcsec) 
or larger (3 arcsec) search radius does not affect the 
statistics of our final sample. Of these 378 sources, 
we only considered the 144 galaxies with a MUSIC redshift 
between 0.7 and 1.5, for which the Ha emission line is ex¬ 
pected to fall in the WFC3-G141 grism spectral range. The 
spectra of these sources have been visually inspected to dis¬ 
card faint or saturated objects or sources without emission 
lines. We excluded also sources that are strongly contam¬ 
inated by nearby objects, carefully evaluating by eye 
the 2D and ID spectra of each galaxy, also con- 
sidering the quantitative contamination model of 


Brammer et al. (2012|. Galaxies whose spectra are 
strongly contaminated were discarded from the fi¬ 
nal sample. To avoid mis-identification of the lines, we 
also discarded objects with \zmusic — zsd-hstI > 0-2. 


GOODS-S field 




21005 30-HST 
591 PACS-MUSIC 
201 PACS-MUSIC outside 3D-HST field 
378 PACS/3D-HST counlerports 
79 This work subsomple 


9755 PACS-MUSIC sources hove o 3D-HST counterport 


53.1 53.2 

R.A. (decimol degrees) 


Fig. 1. Overlay showing the spatial distribution of the 
PACS/MUSIC/3D-HST source sample. 


2.4.1. The X-ray luminous AGN population 


We used the Chandra 4Ms catalog from Xue et al. (2011 1 to 
identify and disc ard the AGN from our sample. Following 
Xue et al. (20111, we classified as AGN only those objects 
with an absorption-corrected rest frame 0.5-9 keV luminos¬ 
ity Lx -ray ^ 3 X 10"^^ crgjs: we found that 8 objects match 
this criterion. We cannot confirm the presence of the AGN 
relying on the 3D-HST data alone, as they do not have 
neither the required resolution to deblend the [Nil] and Ha 
lines, nor the adequate spectral coverage to detect other 
emission lines entering the BPT diagram (Baldwin et al. 
TM]|, e.g. [OIIIlsooT and H/3. 


2.5. The final sample 

To summarize, since we wanted to study galaxies with de¬ 
tections in the far-IR and in the Ha emission line, we con¬ 
sidered only the 144 galaxies in the range 0.7 < z < 1.5 
from the full sample of PACS objects with 3D-HST coun¬ 
terparts. The spectra of these sources were then visually 
inspected to discard the faint, noisy, saturated or strongly 
contaminated ones, removing AGN as explained in the pre¬ 
vious section. 

The final sample includes 79 sources, i.e. ~55% of the 
144 PAGS-selected objects at 0.7 < z < 1.5 with a coun¬ 
terpart in the 3D-HST catalog. The sample galaxies have 
observed Ha luminosities LiHa.obs € [1-4 x 10^^ —5.6 x 10"^^] 
ergjs and bolometric infrared luminosity Lip S [ 1 . 2 x 10 ^° — 
1.3 X 1012] L0. An overlay showing the spatial distri¬ 
bution of our source sample is reported in Figure 

m 

Due to the adopted selection criteria, our sources are 
not fully representat ive of the whole Main Sequence (MS ) 
population at z ^ 1 ( Elbaz et al. 12007 , Noeske et al.j 2007| ) . 
This effect is visible in Figure ^ showing the distribution 
of our sample in the M*-SFR plane with respect to the 3D- 
HST galaxies at z € [0.7 — 1.5]. The stellar masses and SFR 
for the overall 3D-HST po pulation are derived from SED 
fitting (Skelton et al.||2014|| . The SFR for the galaxies ana¬ 
lyzed ill this work are instead estimated from the IR+UV 
luminosities and the stellar m asses M., are computed us¬ 
ing the MAGPHYS software (da Gunha et al. [2008 1 , as 


detailed in Sections |4.1| and |4.2| Fig. [^ offers a qualitative 
comparison of the two selections, and has the only purpose 
of showing the effects of the used selection criteria on the 
global properties of our sample. Our objects occupy the 
upper part of the MS of “normal” star-forming galaxies at 
z ^ 1, due to the adopted far-IR selection (corresp ondng 
to a selection in SFR, cfr. Rodighiero et al. 20141. Only 
one galaxy lies 4x abov e the MS, i.e. the starburst region 
(Rodighiero et al. ]MI| . i lowever, the presence of this out- 
lier does not influence our results, hence we did not remove 
it from the sample. The SFR and mass ranges spanned by 
our sample are 3 < SFR/ij+yy < 232 M 0 /yr, 2.6 x 10® < 
M* < 3.5 X 1 O^^M 0 , respectively. 

Figure [^ shows the dependence of the total infrared lu¬ 
minosity as a function of the redshift for the PAGS-MUSIG 
sources that are located inside the area observed also by 
the 3D-HST survey: 97 % of these objects have a 3D-HST 
counterpart. In Figure]^ the 79 sources of our final sam¬ 
ple, namely the PAGS-MUSIG/3D-HST sources with an Ha 
emission and a “clean” spectrum, are also highlighted. 


3. Spectral analysis 

We extracted the ID spectra by collapsing the 2D spectra 
inside specific apertures, adapting an IDL procedure which 
also applies the calibration files by Brammer et al. (2012). 
Figure [^ shows the grism spectrum and the extracted ID 
spectrum for one of the sources in our sample. The upper 
panel displays the cutout of the grism exposure: the 2D 
spectrum of the object is located in the central part of the 
frame. 

The extraction of the ID spectrum (lower panel of Fig¬ 
ure S is done inside a frame region called “virtual slit” 
(insiae the two red horizontal lines in the upper panel of 
Figured. The width and the position of the virtual slit are 
adjustea ad hoc on each object to maximize the S/N ra¬ 
tio and to minimize the contribution to the flux from other 
spectra. The flux contribution from other sources is negli¬ 
gible in the extracted spectra, as we defined the position 
of the virtual slits to minimize this effect. Furthermore, 
for the majority of galaxies in our final sample, this 
confusion usually affects the flux at all wavelengths, 
so the contamination, if present, is mostly removed 
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Fig. 2. Location of our sample (red filled circles) in the SFR- 
M* space, compared with the distribution of the all 3D-HST 
galaxies (black dots) in the same redshift ra nge. The blue solid 
line represents the Main Sequence at 2 ~ 1 ( Elbaz et al.pOOT l 
and the dot-dashed blue line is the 4xMS. The grey dashed line 
represents the SFR corresponding to the 3a flux limit at 160 
(/a=2.4 mjy) for a typical main-sequenc e galaxy, derived 
using the median SEDs of|Magdis et al. ( 2012||. The inset shows 
the redshift distribution for this work sample (red curve) com¬ 
pared to the complete 3D-HST sample in the same redshift range 
(black curve). The lower panel on the left is the distribution of 
the observed Ha luminosities for our sample (not corrected for 
dust attenuation), the lower panel on the right shows the distri¬ 
bution of L/_r. 


r 


.m, 

V t?/ 


I • PACS-MUSIC somple 

1 • PACS-MUSIC/3D-HST (97%) 

1 • This work subsomple (55% of ze[0-7-1.5]) 

- . 

_ 




01 2 3 4 5 6 7 


Fig. 3. Total infrared luminosity Lm as a function of redshift 
for the PEP sample. The magenta Hlled circles are the 378 PEP 
sources with a counterpart in the 3D-HST observations. The 
blue filled circles highlight the 79 sources of our final sample, 
located in the redshift range 2 € [0.7 — 1.5] 


by the flux correction to the spectrum as discussed 
in Sect. |3.2[ For 9 sources out of 79, we found that 
the flux is strongly contaminated at the edge of the 
spectrum but the spectral range around the Ha and 
the closeby continuum are not affected by this con¬ 
tamination. We accounted for this problem in these 
9 sources by carefully selecting the portion of the 
continuum near the Ha emission when rescaling the 
spectrum to the observed broad-band photometry 
(i.e. to derive the aperture correction factor) and 


excluding the part of the continuum affected by the 
contaminat ion. 




1.1x10^ 1.2x10^ 1.3x10* 1.4x10* 1.5x10* 1.6x10* 1.7x10* 

Observed wovelength [A] 

Fig. 4. Upper panel: image of the 3D-HST bi-dimensional spec¬ 
trum for source 4491 (ID-MUSIC), as displayed by the IDL 
routine. The red horizontal lines dehne the “virtual slit”. Lower 
panel: ID spectrum of source 4491, obtained collapsing the 2D 
spectrum along the columns inside the virtual slit. The flux 
is in erg!secjcrn^jK. The black thicker histogram is the ob¬ 
served spectrum, the red curve is a spline interpolation to the 
observed spectrum. The grey thin histogram, almost over¬ 
lapping the black one, represents the "decontaminated 
spectrum", i.e. the contamination model by |Brammer| 
et al.|| ||2012|) subtracted to the observed spectrum. The 
overlap between the observed and the decontaminated 
spectrum highlight the fact that the contamination flux 
in the region defined by the virtual slit is negligible. 
The vertical blue line highlights the position of the Ha emission 
while the pale blue lines correspond to the position of the [Nil] 
doublet. 


3.1. Ha fluxes and redshift measurements 


We measured the redshift and the integrated Ha flux from 
the extracted ID spectra by fitting the emission lines with 
a Gaussian profile based on the IRAF tool splat, even if the 
spectral line profile is dominated by the object shapes (the 


so-called “morphology” broadening, Schmidt et al. 20131 


and the grism line shape may not be well described by the 
Gaussian profile. For part of the sample sources we 
re-measured the fluxes by simply integrating the 
area under the line, without fitting any profile, and 
found that the flux measurements are in agreement 
with those obtained with the Gaussian fit. The errors 
on the redshift and on the Ha integrated flux were esti¬ 
mated using Monte Garlo simulations with random Gaus¬ 
sian noise. The measurements of the observed Ha luminosi¬ 
ties and the redshifts are listed in Tab. |B.l| in Appendix [B| 
The redshifts are distributed between 2 ~ 0.7 and 1.5 as 
shown in Figure as a result of the combined transmis¬ 
sion of the grism G141 + F140W Alter. Our measure¬ 
ments are in good agreement with the pre-existing red¬ 
shift measurements. Figure shows the excellent correla¬ 
tion between the MUSIG spectroscopic (black dots) and 
photometric (red dots) redshifts and our redshifts mea¬ 
sured from 3D-HST spectra. The median absolute scat¬ 
ter Az = \z 3 d-hst — zmusicI is 0.0040, while the me- 
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dian relative scatter /Az/{1 + z^d-hst) — 0.0022. We also 
compa re our redshif t s wit h those from the recent cata¬ 
log of Morris et al.| (20151, measured also from 3D-HST 
spectra: also in this case a very good agreement is found 
(Az = \zsD-HST — ZMorris\ — 0.0051, Az/(1-f Z3 D_//St) — 
0.0027). 



Fig. 5. Redshift distribution for the 3D-HST sources in the 
GOODS-South field. The sources are distributed in a redshift 
interval z G [0.65 — 1.53], according to the features of the 
WFC3/G141 grism. 



0.6 0.8 1.0 1.2 1.4 1.6 

^30-HST 


Fig. 6. Upper panel: relation between 3D-HST spectroscopic 
redshifts (i-axis) and MUSIG redshifts. In the inset, the distri¬ 
bution of the absolute scatter Az={z3d-hst-zmusic) is shown. 
The standard deviation for this distribution is (t= 0.031. The 
red histogram highlights the absolute scatter for photomet¬ 
ric MUSIG redshifts. The blue histogram represents the abso¬ 
lute scatter between our measurements and the redshifts from 


Morris et al. (2015 I. Lower panel: relative scatter {z3d-hst- 
zmusic) / (I+^sd-hst). The data points in red are the photo¬ 
metric redshifts in the MUSIG catalog, the black dots are the 
sources with spectroscopic redshift also in the MUSIG catalog 
from ground based measurements. The blue open diamonds are 
the redshifts measured by Morris et al.. 


3.2. Flux correction by SED scaling 


The position and the width of the virtual slit used to ex¬ 
tract the ID spectrum were adjusted for each object in 
order to maximize the S/N ratio and minimize the contam¬ 
ination from nearby sources. We applied corrections to the 
measured Ha fluxes to account for flux losses outside the 
slit. For each ID spectrum the correction factor was de¬ 
rived as the ratio between the average continuum flux in 
the ID s pect rum (Fgpec) and on the galaxy SED (Fsed, 
see Sect. 4.1 for details about the SED fitting), measured 
in the same wavelength range. The average value of the 
correction factor is FsEo/Fspec ~ 1.41, ranging from 


FsEo/Fspec ~ 0.39 up to FSED/Fspec 8.24 . 

Once this flux correction is applied, all the objects show 
an excellent agreement between the ID spectrum and the 
near-IR part of the SED. An example is shown in Figure 



Fig. 7. Spectral Energy Distribution from MAGPHYS (black 
curve) with the observed photometry (blue filled circles) and the 
3D-HST spectrum (red line). The 3D-HST spectrum is corrected 
for the flux loss outside the slit, as described in the text. The 
plot also reports the MUSIG ID of the source and its redshift, 
measured from the near-IR spectrum. The inset shows a zoom 
on the 3D-HST spectrum, close to the Ha emission line. 


4. Derivation of the main physical quantities of the 
sample galaxies 

4.1. SED-fitting 


We fitted the Spectral Energy Distributions with the MAG¬ 
PHYS package ( Multi-wavelength Analy sis of Galaxy Phys¬ 
ical Properties, da Gunha et al. 2008). This software has 
the main advantage of fitting the whole SED from the UV to 
the far-IR, relating the optical and IR libraries in a phys¬ 
ically consistent way. To compute the SEDs, we ran 
MAGPHYS in the default mode, using the stellar- 
population synthesis models of|Bruzual & Chariot] 
( 2003 ). For the SED fitting we used 21 photometric bands, 
i e. from the near-UV (GALEX) to the SPIRE 500 /rm. The 
Alters used for the SED fitting are listed in Table The 
errors of the photometric data points were set to be 
10% of the measured flux and were then adjusted 
ad hoc. The best-fit SEDs were obtained by fixing the red¬ 
shifts to the spectroscopic values that we measured from the 
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3D-HST spectra. Figure shows an example of a best-fit 
SED as output by MAGPHYS. 

In addition to the best-fit SED, MAGPHYS returns sev¬ 
eral physical parameters of the observed galaxy together 
with their marginalized likelihood distributions. We used 
in particular the stellar masses M* and the bolometric in¬ 
frared luminosity L/^j, considering the values at the 50*^ 
percentile of the likelihood distribution. The errorbars 
at 68% for these quantities were derived from the 
PDFs c omputed by M AGPHYS. Since MAGPHYS 
adopts a Ghabrier (2003) IMF, we converted the stellar 
masses to a Salp eter IMF by multi plying by a constant 
factor of 1.7 (e.g Cimatti et al.|[2008 |. 


Table 1. Filters used for the SED fitting and respective effective 
wavelengths Ae// (mui)- 


Filter 

Ae// [/rm] 

GALEX NUV 

0.2310 

U 

0.346 

AGSF435W 

0.4297 

AGSF606W 

0.5907 

AGSF775W 

0.7774 

AGSF850LP 

0.9445 

ISAAGJ 

1.2 

ISAAGH 

1.6 

ISAAGK 

2.2 

IRAGl 

3.55 

IRAG2 

4.493 

IRAG3 

5.731 

IRAGl 

7.872 

IRS 16 

16 

MIPS24 

23.68 

PAGS70 

70 

PAGSlOO 

100 

PAGS160 

160 

SPIRE250 

250 

SPIRE350 

350 

SPIRE500 

500 


4.2. Measuring the SFRs 

We derived for all the sources in the sample the SFR from 
several indicators, thanks to the availability of a wide multi¬ 
wavelength photometric coverage and the near-IR spectra. 

In particular, we computed the SFRs from the Ha 
luminosity, the UV luminosity and the IR luminos- 
i ty, ad opting the classical calibrations of |Kennicutt| 
(]1998|). 

The Ha luminosity is derived from 3D-HST near-IR 
spectra. The meas ured Ha flux is aperture corrected as de¬ 
scribed in section 3.2 while the [Nil] contribution was re- 
mo ved by scaling down the Ha flux by a factor 1.2, follow¬ 
ing Wuyts et al. (2013). The contamination by [Nil] could 
produce some uncertainties in the estimate of I^Ha since 
the ratio [NH[/Ha may vary object by object as a func - 
tion of the stellar mass and the redshift (Zahid et al. ]M4| . 
To validate our choice of using a constant factor, we 
then computed the [Nil] correction factor as a func¬ 
tion of M*, using the evolutionary parametrization 



Fig. 8. Example of a Spectral Energy Distribution in output 
from MAGPHYS. The black solid line is the best-fit model to 
the observed SED (data-points in red). The blue solid line shows 
the unattenuated stellar population spectrum. The bottom panel 
shows the residuals (LfT 


of the ISM metallicity by Zahid et al. fl2013). We 
found that the use of a more accurate iNiTr correc¬ 
tion does not change our results about the dust ex¬ 


tinction corrections (i.e. Figures |10| and 11). Results 
of our later analysis (like the /-factor and the com¬ 
parison between SFR indicators) are not affected 
by the choice of a mass-dependent [Nil] correction 
or of a constant scaling factor. 

In order to account for both the unattenuated and the 
obscured Star Formation, we measured the “total” SFR 
combining the L/^ an d the observed Lifiop (the same pro¬ 
cedure was adopted by Nordon et al.|]20l^ : 

SFRir+uv = (1.4Li6oo + 1-7Lir) x 10-^°[Lq] . (1) 

We estimated the luminosity at \rest-frame = 1600 A from 
the best-fit SED, and used the IR luminosity L/jj derived 
by MAGPHYS (cfr. Section |^. The SFRir+uv is the 
best estimate of the SFR for our objects, considering the 
wide photometric coverage in the far-IR, and will be used 
as a benchmark to verify the validity of our dust extinction 
correction. 

The errors for the various estimates of the SFR are com¬ 
puted via MG simulation by randomly varying Fua, Ligoo 
and Fir within the errorbars, assuming that these quanti¬ 
ties are Gaussian distributed. 


5. Dust extinction corrections 

In order to precisely evaluate the SFR and to reconcile the 
various SFR estimates, we need to accurately quantify the 
effect of dust on the integrated emission of our galaxies. 

According to Calzetti (2001), the action of dust on 
starlight for starburst galaxies in the local Universe 
can be parametrized as: 


FobsW = F„t(A)xl0-°'^^" = 

( 2 ) 

where Fobs (A) and Fi„t(A) are the dust-obscured 
and the intrinsic stellar continuum flux densities. 
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respectively. A\ and 'EgtariB — V) are the dust atten¬ 
uation and the color excess on the stellar continuum 


respe c tively , while fc(A) parametrizes the Calzetti 
et al. (1994) starburst reddening curve. 


!Ss alrea dy stated in the introduction, |Calzetti| 


et al. (20001 found that exists a differential attenuation 
between the stellar continuum and the nebular lines by 
measuring both the hydrogen line ratios and the contin¬ 
uum reddening in a sample of local star-forming galaxies. 
This differential attenuation is parametrized by a /-factor 
of 0.44, as defined by : 


E^tariB - V) = f X Er,eb{B - V) 
with / = 0.44 ± 0.03 . 


( 3 ) 


We note that in the original calibration of eq. (3) two 
different reddening curves were use d to measure the con- 
tinuum and the nebular extinction (Fitzpatrick 1999| and 
Calzetti et al. 2000[ respectively). If the Calzetti redden- 
ing curve is used to measure both the nebular and contin¬ 
uum extinction components, as done in this work, the local 
/-factor becomes /=0.58, instead of the canonical /=0.44 
( Pannella et al.|[2014 Steidel et al. |2014 1. Hence, hereafter 
we refer to /=0.58 as the “local” /-factor. 

Equation 3 implies that the ionized gas is about two 
times more extincted than the stars. The applicability of 
such relation in the high redshift Universe is still not clear, 
as already mentioned in Sect. 1. Several authors (e.g. Pan¬ 


nella et ai. 2014[|Kashino et al.|2013[|Erb et al.|2006[ Reddy 


et al.|20l0 Price et al.|2014 Wuyts et al.|2013[ | found that 

the use of this equation implies to overestimate the intrin¬ 
sic line flux, i.e. the Ha line is less attenuated in the high 
redshift galaxies with respect to the local Universe. 

To shed some light on this issue, we investigated the dif¬ 
ferential extinction between nebular lines and continuum 
at high-z by using the ratio between the observed SFRijc 
and SFRf/y, thus uncorrected for dust attenuation. To val¬ 
idate the derived dust correction we compared our measure¬ 
ments of SFR^fa to the SFR/jj+c/r^, that we assume as the 
“true” SFR estimator. 

In the following section we present the formalism and 
the methods to derive the Ha differential attenuation. 


5.1. Differential extinction from Ha to UV-based SFR 
indicators 

We use the Ha and UV luminosities, both uncorrected for 
dust extinction, to derive the extra correction factor associ¬ 
ated with the nebular lines as a function of the color excess 
on the continuum emission. For the reader’s convenience, 
we report in this Sect, the relevant formalism. 

We assume that the total Star Formation Rate SFR^oi 
is proportional to luminosity: 

SFRtot = const X Li„t(A) (4) 

where, in the case of UV-optical wavelengths, the intrin¬ 
sic luminosity Li„t(A) is related to the observed Lo{,s(A) 
through the attenuation Ax (eq. [2 1. Combining equations 
and 1^ we obtain for the observM SFR: 

SFRuncorr = COUSt X io6s(A) = X SFRtot ■ (5) 

Considering the UV and the Ha SFR tracers and 
assuming that the intrinsic SFRs in absence of dust 
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attenuation are the same (SFR//q, = SFRf/y), we 
obtain that the ratio between the observed SFRs 
is related to the differential extinction in the Ha 
emission line: 


SFRna.uncorr X SFRtot 

SFRuv,uncorr x SFRtot 

^ ^Q-0.4E,t^r.(B-V)x[!Ft^-k{UV)] 


or in logarithmic units 

, ,SFRh a.uncorr-i a 

log[ -1 = —0.4 


SFRifv, 


uncorr 


k{Ha) 

f 


k{UV) 


( 6 ) 

X Estar{B - V) . 

( 7 ) 


We can then quantify the differential extinction 
of nebular lines through a linear fit in the plane 


from eq. j? 


E,tar{B-V),log[§^^^ 

We stress here that in this analysts we impose 
that the intrinsic Ha and UV luminosities are trac¬ 
ing the same stellar populations and this condition 
could be in principle not satisfied, as these two 
SF tracers are sensitive to different star formation 
timescales. In fact, the UV luminosity is dominated 
by stars younger that 10 ® yr while only stars with 
masses higher than 10 Mq and lifetimes lower than 
~ 20 Myr contribute significantly to the ioniz ation 
of the HII r egions, i.e. to the Ha luminosity (Ken- 
tt||1998 ). However, our assumption seems to be 


reasonable from a statistical point of view, as we 
would need to observe a galaxy when the formation 
of stars with ages < 10^ yrs (i.e. the Star Formation 
traced by L(Ha)) is already exhausted, while the 
formation of stars with ages 10^ < t < 10® yr is al¬ 
ready in place (i.e. sources that produce L(UV) but 
do not contribute significantly to L(Ha)), to have 
L(Ha) substantially different from L(UV)). 


5.2. Measurements of the continuum extinction 


After we have parametrized the /-factor and measured the 
ratio SFR//Q,_„„corr/SFR( 7 v^„„corr, we need to estimate the 
continuum color excess Estar(B-V). 

Taking advantage of the far-IR measurements available 
for our sample, we derived the continuum atten uation Aipx 
from the ratio L/^ {Ijuv (Nordon et al.|[2M2|: 


Airx = ‘l.blogl 


SFRir+uv , 
SFRijv 


( 8 ) 


We can independently infer the continuum attenuation 
also from the UV spectral slope /3 because it is a sensi¬ 
tive indicator of dust attenuation. The intrinsic shape 
of the UV continuum spectrum for a star-forming 
galaxy is nearly flat in F(A), even considering two 
extreme situations: an instantaneous burst of star 
formation and a constant star formation rate. For 
the case of the instantaneous burst, the usually as¬ 
sumed burst duration is typically ~ 20 Myr. During 
the first 2x10^ yr, the stars contributing to the flux 
in the range A G [1200 — 3200] A have had no time to 
evolve off the main sequence, so that the shape of 
the intrinsic spectrum in the UV wavelength range 
of interest has not changed. Considering instead a 
region of constant star formation, the continuous 
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generation of new stars keeps the shape of the UV 
spectrum roughly constant. In conclusion, we can 
reasonably assume that the intrinsic shape of the 
UV spectra of star-forming galaxies is constant so 
each dev iation from this intri nsic shape is produced 
by dust (Calzetti et al. 1994). 


We then computed the continuum atten uation at 1600 
A (Ai eoo) from /3 using the calibration from Meurer et al. 
(19991, derived on a local sample of starburst galaxies: 


^1600 — 4.43 + 1.99^ . 


(9) 


Here the UV-slope parameter /3 is defined as the linear in¬ 
terpolation of the observed spectrum (or in the lack of it 
from the photometric data) in the rest-frame wavelength 
range A € [1250 — 2600]A. The applicability of the Meurer 


strated by several authors (e.g. Reddy et al. 

2012 

Buat 

et al. 

2012 

Talia et al. 

2015D. 


In our case, since we lack observational data in this 
wavelength interval for the majority of the sample, we es¬ 
timated /3 from a linear interpolation of our best-fit MAG- 
PHYS spectral mo del discussed in Sect. 4.1, as done e.g. in 
Oteo et al. (20141. To test the robustness of this slightly 
model-dependent estimate of the UV-slope, we also com¬ 
puted /? from the UV rest-frame photometry, when available 
(i.e. for 13 objects), and found that the two measurements 
are in good agreement. For more details see Appendix [Aj 
The attenuation is related to the color excess E s tarfB- 
V). Assuming the reddening curve of Calzetti et al. ]( [2000t 
we have: 


Estar,IRx{B — V) 


Airx 

k{UV = 1600) 


( 10 ) 


Estar,p{B - V) 


^1600 

k{UV = 1600) ■ 


( 11 ) 


1996|, of 

total scatter 


These two quantities are consistent to each other, as dis¬ 
played in Figure]^ The correlation between Es(ar,/ 3 (B-V) 
and Estar./fljf (B-V) is very good, confirming the reliabil¬ 
ity of our estimate of continuum dust extinction, and has 
an intrinsic scatter, i. e. not accounted for by the ex¬ 
perimental errors (Akritas &: Bershai 
o’mtr- ~ 0.061, that represent 95% of the 
(ctot ~ 0.064). Estimates of amtr a nd atot are obtaine d 
with the IDL routine mpfit.pro (Markwardt 2009). 

Note that, since our UV-slope /3 is measured from SEDs 
computed with MAGPHYS, which consistently accounts for 
both the UV and far-IR emission, this could partially ex¬ 
plain the high correlation degree between Estaj.,/ 3 (B-V) and 
EstaT-,//jx(B-V). In any case, this comparison is useful to 
verify the accuracy of our continuum extinction measure. 

Both measurements of Estaj.(B-V) will be considered in 
the subsequent analysis to evaluate the differential extinc¬ 
tion, i.e. the /-factor. 


The errors for Estar,/ 3 (B-V) were derived by propagating 
those on /3 (see Appendixfflfor the calculations of the errors 
on the UV slope /3), while the errorbars for EstarjRx{B—V) 
are estimated via MG simulations, by varying L/^j and Lj/y 
within the la errorbars. 



[mag] 


Fig. 9. Comparison between the continuum extinction 
Eatar{B — V) derived from the UV spectral slope /3 (x-axis) 
and the continuum extinction derived from the infrared excess 
IRX (y-axis). The red line is the 1:1 correlation relationship. In 
the lower part of the plot on the right are indicated the typical 
size of the errorbars. 


5.3. Extra extinction from the f-factor 

Following the formalism of section |5.1| we evaluated the 
/-factor as a function of the color excess on the contin¬ 
uum by using the ratio SFRHa,uncorr/SFRuv,uncorr and 
the two measurements of the color excess on stellar con¬ 
tinuum, Estar./3(B-V) and Estar,/itx(B-V). To derive the 
/-factor we fitted a linear function y = s x x to the data 
points through a minimization of a weighted (i.e. a 
least square fit obtained by weighting each data 
point for its error, both in the x and y axes), where 
y ^^9(.^^Brra,uncorr/^ERijyy^^(.Qj.j.^ and X —Egia,.(B-V). 
To obtain this fit we used the routine mpfit.pro, and took 
care of checking the result by a MG simulation by insert¬ 
ing a scatter in the x-coordinate value for each datapoint. 
The /-factor is then finally obtained from the slope s of the 
best-fit line as 

0.4 X k{Ha) 

^ 0.4 X k{UV) -s' ^ ^ 


The uncertainty for the /-factor was estimated by propa¬ 
gating the error on the best-fit parameter s. 

Figure 10 reports SFR//Q,«ncorr/SFR[/y_ 

U 


as a 


function of the continuum color excess. The red solid line in 


Figure 10 corresponds to Equationl^with /=0.93 (± 0.065) 
that is the best-fit value from our data distribution. Hence 
we found that the /-factor required to match SFR//q, and 
SFR[/y in our sample is larger than the local value (/=0.58, 
blue dashed line in Fig. [Io| an d turns out to b e also h igher 
than the values computed by Kashino et al.j (2013) on a 
sample with slightly higher z and different selection crite¬ 
ria (/=0.69 and 0.83, black dot-dashed and green dot-dot- 
dashed lines in Fig. [To] respectively). 

The value of the /-factor does not change using either 
Estar -./3 (B-V) or Egtar./flx (B-V): this result is in line with 
the high correlation of these two measurements, as seen in 

Fig-i 
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Fig. 10. Ratio of Ha to UV-based SFRs (not corrected 
for dust extinction) as a function of Estar{B — V) derived 
from the IRX ratio . The lines are the e quati on [7| with dif- 
faent values of / from Kashino et al. (20131 and Calzetti et al.| 
(20001, as the legend indicates. The grey shaded area marks the 
conhaence interval for our estimate of the /-factors. 


5.4. Testing the dust correction: comparison between SFR^a 
and SFRir+uv 


We tested the reliability of our extinction correction in 
Figure 11 where we compare the SFR/^j+f/y with the 


SFRjjq,, corrected for dust extinction using both the pre¬ 
scription derived in this work (i.e. EistarjRX (B- V) and 
/=0.93, upper panel), and that derived by jCalzetti 
et al.| (|2000|) for local galaxies (Estar,/itx(B-v) and 
/=0.58, lower panel). The results show that a bet¬ 
ter agreement between SFR//q, and SFR/i{+[/v is ob¬ 
tained by applying our correction factor, with a me¬ 
dian ratio SFR^fQ/SFR7;j+[/y=0.88 and a median rel¬ 
ative scatter |SFR/77y[/y-SFR7fa|/(l+SFR/77y(7y)=0.38. 
On the contrary, the Calzetti et al. prescription would lead 
to overestimate the SFR^/a by a factor up to 3.3 above 
SFR/i7+c/y >50 M0/yr. We emphasize that our recipe 
for the nebular dust attenuation seems to “fail” for the 
tail of objects with SFRjp+uv F [10 — 2OM0/yr and 
M* ~ 1O^°M0 (light blue points in Fig. JH). For these 
sources our dust correction underestimates the SFR/^q, with 
respect to SFR/Tj+j/y, while using the Calzetti prescription 
the two SFRs are in slightly better agreement. This could 
imply that the /-factor is an increasing function of SFR and 
M*. 

However, the disagreement between SFR indicators at 
SFR/7j_|_77y < 20 M0 /yr can be due also to an overestimate 
of SFR/7j_|_7/y rather than a problem of the dust correction. 
At lower M* and SFRs in fact, the contamination to the 
heating of dust by an ol der stellar population ( the so-called 
“cirrus” component, e.g. Kennicutt et al. ]200^ became not 
negligible and can cause to overestimate L/tj. Finally, the 
trend seen in Figure E] could also be a consequence of 
the selection criterion used in this work, considering that 
SFR/7j+[/y ~ 15M0/yr is the value that corresponds to 
the 3(7 flux limit at 160 pm at z ^ 1 a s estimated using 
the median SEDs of Magdis et al. (2012|| (see also Fig.|2|of 
section 2.5). 
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5.5. Ah a '/s M, 


Figure shows the relationship between the attenuation 
Ana and stellar masses. Here we converted the continuum 
attenuation E star.iRxjB — V) to Ah o by using f=0.93 and 
the k(A) from Calzetti et al. (2000). Despite the high dis¬ 
persion of the data, we can observe an excess of Ana at 
the highest masses w ith respect to the local relationship of 
Garn & Best| (|2010l (green dashed line in Figure [T^. We 
report also the relationship of|Kashino et al. (2013|), oerived 
at z ~ 1.6 (blue solid line). The grey tilled dots of Figure 12 
represent the median values of Ana obtained in three M7 
bins as reported in the figure. The errorbars are the median 
scatter of Ana and M* in each bin. Within the large scat¬ 
ter in the data, the median values in each bin are consistent 
with both the local and the higher-z relat ionship, in agree¬ 
ment with the results of Ibar et al. ( 2013[ ). Of course, there 
may be selection effects in operation in the graph disfavor¬ 
ing the detection of galaxies with high values of Anai that 
may be more important for the lower mass objects with 
intrinsically faint Ha flux. Our combined selection requir¬ 
ing detection of the Ha line may be somewhat exposed to 
such an effect. To confirm if an evolution with z ex¬ 
ists in the dust properties of star forming galaxies 
it would be necessary to have deeper observations 
for the most attenuated galaxies. This aspect of the 
analysis will be investigated in a further paper. 



Fig. 12. Ana as a function of stellar masses. The black filled 
circles are the sources at z < 1 while red filled squares are objects 
with z > 1. The attenuation is computed from EstarjRX (B-V) 
by using f=0.93. The grey filled circles are the median values 
of Ana in three mass bins (M* < 3 x 10^° Mq, 3 x 10^° < M* 
< 1.7 X 10^^ M©, M* > 1.7 X 10^^ Mq) and the errorbars are 
the median scatter for Aho and M* in each mass bin. The two 
vertical dashed lines highlight the mass bins (M* < 3 x 10^° Mq, 
3 X 10^“ < M* < 1.7 X 10“ Mq, M* > 1.7 x 10“ Mq). 


5.6. Main sequence at z ^ 1 

In Figure |13| we report an updated version of the 
relation between the stellar mass and the SFR, al¬ 
ready presented in Figure where SFR is now 
computed from the Ha luminosities, corrected for 
dust extinction according to our reference dust cor¬ 
rection. As already mentioned in Section |2.5[ our 
sources do not span the overall range covered by the 
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star-forming Main Sequence population at z ^ 1, 
since our far-IR selection favor the detection of 
sources with SFR>10-20 Mq/yc. However, we note 
that above stellar masses on the order 3 x lOi°M0, 
our sample is basically representative of a mass- 
selected sample that traces the underlying MS at 
the same redshift. 



Fig. 13. SFRrfa versus M* with hna corrected by dust with 
our recipe. The blue solid line reports the MS relationship from 


Elbaz et al. (20071 at z = 1 and the blue dashed line is the 
corresponding Ax MS. The black filled circles are the sources 
at z < 1 while the red filled squares are the objects at z > 1. 
The median errors for M* and SFRrfa are reported in the lower 
part of the plot on the right. 


6. Caveats: the many faces of the f-factor 

Values of the /-factor estimated in the literature range from 
0.44 to ^ 1 , as shown in Table and in this section 
we briefly discuss the likely origins of these differences. An 
/-factor less than unity is currently interpreted as implying 
more average obscuration affecting the line emitting regions 
(HIT regions) compared to the average obscuration affecting 
the hot stars emitting the UV continuum, i.e., / yf 1 would 
be the result of a geometrical effect with the spacial dis¬ 
tribution of line emitting regions being different from that 
of continuum emitting stars. Actually, the reality may be 
somewhat more complicated. 

First of all, the /-factor can be derived in two radically 
different ways: either by comparing two extinctions or two 
SFRs. In the former case the Ha extinction derived from the 
Balmer decrement is compared to the extinction in the UV 
as derived from either the UV slope or from Equation (10). 
Then one has to adopt one specific reddening law k{\). 
In the latter case, the /-factor is estimated by enforcing 
equality between the SFR derived from the Ha flux with 
the SFR derived from another indicator, such as SFR(UV), 
SFR(UV+IR) or SFR from SED fitting, etc. 

When the /-factor is derived by comparing Ha and UV 
extinction, then the result depends on the adopted red¬ 
dening law fc(A), hence it reflects at once both the men¬ 
tioned geometrical effect and possible departures from the 
adopted reddening law. As well known, the reddening law 


is not universal, not even within the Local Group. These 
two contributions to determine the value of / can barely 
be disentangled. When the /-factor is derived by forcing 
agreement between SFR(Ha) and the SFR from another 
indicator, then it becomes a fudge factor to compensate for 
the relative biases of the two SFR estimators, neither of 
which will be perfect. 

Moreover, the measured Ha flux is subject to dust ab¬ 
sorption in two distinct ways: first, each Ha photon has 
a probability escape the galaxy, but, second, 

dust absorption in the Lyman continuum reduces the num¬ 
ber of produced Ha photons by a factor , 

where ALyman-cont. is the extinc tion in the Lyman con¬ 
tinuum. With few exceptions (e.g. Boselli et al.||2009 ), this 
second aspect is generally ignored, in the hope that an em¬ 
pirical calibration of SFRs may subsume in it also this part 
of the involved physics. In any event, the resulting /-factor 
depends on the actual extinction law of each galaxy, extend¬ 
ing from the optical all the way to the Lyman continuum, 
on the geometry of the emitting regions, as well as on the 
relative systematic biases in the relations connecting SFRs 
to observables. 

Besides these aspects, derived / values can also depend 
on the specific galaxy sample from which it is derived. 
For example, in our approach the far-IR selection crite¬ 
rion leads to select objects with strong levels of dust ob¬ 
scuration in the UV. The other requirement is to select ob¬ 
jects with a strong Ha emission {Fohs{Ha) > 2.87 x 10“^^ 
ergjs and S/N higher than ~ 3), hence with lower lev¬ 
els of Ha attenuation. The two selection criteria partially 
conflict with each other and combined favor sources with 
E„eb(B-V) ~ Estai.(B-V). To understand how the selection 
criterion influenc es the results we can compare our analysis 
with the work of Kashino et al. (20131. Our analysis was 
performed f ollowing the same approach of |Kashino| 
et al. (|2013) but our sample has different selection 
criterion and size (168 sBzK galaxies for Kashino, 
79 far-IR sources in this work) and we obtained an 
/=0.93, that is 35% larger than the Kashino result. 
In the case of re st-frame UV-selected galaxies, such as in 
Erb et ah] (2006 1 , this bias favors galaxies with low extinc- 
tion whicn may result in a different / value compared to 
the case of samples including also highly reddened galaxies. 

All of these considerations imply that different indica¬ 
tors lead to different estimates of the /-factor. For example, 
if we consider the ratio SFR^/SFR/i^+c/y, we get (see e.g. 
eq.jT]) 


SFRhc. ,'U 


SFR 


= 10 


— 0.4AHa _ 10 ~*^ ^ 


f 


k(Ha) 


IR+UV 


(13) 


Our data would imply in this case /~ 0.85, however not 
very significantly different from our best-guess of /~ 0.93. 

Also the estimate of F,star{B — V) and its errors influ¬ 
ences the estimate of the /-factor, leading to results that 
can vary from / ^ 0.4 to a value larger than 1. 

The last point to consider is again the reddening law 
fc(A) assumed in the analysis. In the literature one finds very 
different tr ends for A:(A), suc h as the presence of a bump at 
^ 2200A (Fitzpatrick 1^991 for the LMC) or a smoother 
trend (Calzetti et al. |2000| for a starburst galaxy), and 
also very different values for its nor malization Ry = 
Ay/E star(B-V) , that vary from 3.1 (Cardelli et al. 
1989) to 4.05 ( [Calzetti et al. 2000). At redshift ~ 2 
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both galax ies with and wit hout the 2200 A bump appear 
to coexist ( Noll et ^|2009 ). The shape and the normaliza¬ 
tion of the assumed fc(A) strongly influences the value of /, 
both for direct or indirect measurements. As an exercise we 
derived the /-factor using different k{X) expressions in eq. 

Table H summarizes the results, showing that the value 
of the /-factors ranges from ~ 0.7 to ~ 1.2. 

In summary, the /-factor may offer a fair mea¬ 
sure of the relative extinction of emission lines and 
the stellar continuum, still however relying on not- 
completely well defined and understood ingredi¬ 
ents. 


Table 2. Summary of the values obtained for the /-factor as a 
function of the assumed reddening curve k{X): the first column is 
the fc(A) assumed to compute the UV attenuation at 1600 A, the 
second column is the k{X) used to compute the Ha attenuation 
while the last column is the /-factor, obtained from eq. |12[ 


k(UV) 


k(Ha 


/-factor 

Galzetti et al. 

(2000 

1 ( 

Dalzetti et al 

1(2000 

1 0.93 


Ualzetti et al. 

(2000 

1 

bitzpatrick 

|iyyy| 

0.66 


Reddy et al. 

2irmr 

Reddy et al. 

IpTB 

0.96 


Reddy et al. 

2015) 

Fitzpatrick] ([lyyyp 

1.19 


7. Summary and conclusions 


In this work we analyzed the near-IR spectra of 79 star 
forming galaxies at z G [0.7 — 1.5], acquired from the 3D- 
HST survey . The sources were selected in the far-IR from 
the Herschel/PACS observations: the PACS catalogs were 
associated with the 3D-HST observations using the IRAC 
positions of the PACS sources. From the near-IR spectra 
we measured the Ha fluxes and the spectroscopic redshifts 
of the whole sample. 

We computed the SEDs with the MAGPHYS software, 
using data from near-UV to far-IR including the GALEX- 
NUV, the GOODS-MUSIC optical to mid-IR catalog, the 
IRS-16 fim and the far-IR photometry from Herschel PACS 
and SPIRE (i.e. at 70, 100,160, 250, 350 and 500 pm). From 
the SEDs we derived the stellar masses, M*, the bolometric 
infrared luminosities, Pjr, the UV luminosities, CuV: and 
the UV slope, /3. 

We then evaluated the color excess Es(ar.(B-V) from the 
IRX=L/ii;/L(7y ratio and from the UV slope f3 and found 
that these two quantities are in good agreement. In our 
sample the color excess on the stellar continuum ranges 
from Es(ar-(B-V) ~ 0.1 mag to Estar(B-V) ^1.1 mag. 

We computed the dust attenuation on the Ha emis¬ 
sion E„eb(B-V) as a function of Est£,j.(B-V) by comparing 
the SFR/zq, and the SFRj/y, both uncorrected for extinc¬ 
tion. We obtained that the /-factor, which parametrizes the 
differential extinction on the nebular lines, is f=PistariP- 
V)/E„ef,(B-V)=0.93 ± 0.06. This result is consistent within 
the errorbars with the analysis of [Kashino et al. ( 2013| ) 
from the Balmer Decrement and of|Pannella et al. (zOl'fl), 


performed in a similar redshi ft range. Our ana lysis is also 
consistent w ith the results of Erb et al. (20061 and Reddy 
et al. ( 20 To| performed at higher z, as summarized in Table 
^ which collect a list of results from others works. The good 
agreement found in our sample between the SFR/jj+f/y and 
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the SFR//Q, corrected for extinction using our recipe further 
confirm our results. 

From our dust correction we then computed the atten¬ 
uation Ana as a function of M*. We found that Ana is 
increasing with M* and this trend seems to diverge from 
the local relationship: our sources shows an excess of Apa 


with respect to the relationship ofjGarn & Best (|2010 1 for 
M* > 10^^ Mq, suggesting an evolution in the dust prop¬ 


erties of star-forming galaxies with z. 

In conclusion we found that the level of differential ex¬ 
tinction required to match the SFRjjq, with the SFR/jj+j/y 
is lower than in the local Universe, thus Aua ~ Auv for the 
sources in our sample. The value of the /-factor seems to be 
related to the physical properties of the sample rather than 
be dependent on z. The trends of Figures and ^ sug¬ 
gest that the Ha extinction (thus the /-factor) is a function 
of SFR and M*. In particular we notice that in Figure [m 
our dust correction underestimate SFRijc with respect to 
SFR/fl;_|_[/y for sources with SFR < 20 MQ/yr: these sources 
require a lower value of / similar to the local /=0.58. This 
trend could be explained by the two component s model of 
dust sketched in Figure 5 of Price et al. (2014). A galaxy 
with high sSFR is supposed to have an high number of 
OB stars, which are located inside the optically thick birth 
cloud: in this case these massive stars dominate both the 
UV-continuum and the Ha emissions, so the level of atten¬ 
uation for the continuum and the nebular emission will be 
similar {Apy ^ Apa)- On the other and, for a galaxy with 
low sSFR the number of OB stars will be lower, so in this 
case the optical-UV continuum is mainly produced by the 
less massive stars that are located both in the birth cloud 
and in the diffuse ISM: in this case Ana > Apy since the 
Ha emission is produced in a different and more dust-dense 
region with respect to the continuum. 

In this paper we do not examine in depth the conse¬ 
quences of this modellistic approach for the reasons dis¬ 
cussed in the previous section. We defer to a future paper a 
more detailed analysis of the extinction properties of star¬ 
forming galaxies based on the “Intensive Program” (S12B- 
045, PI J. Silverman) with the FMOS spectrograph at the 


Subaru Telescope in the GOSMOS held (Silverman et al. 


2014). 
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Appendix A: Computation of the UV spectral slope 

Appendix A.l: jd from the best-fit SED 

The UV slope is derived from a linear fit to the model best- 
fit SED in the wavelength range A S [1250—2600] A, lacking 
observations in this rest-frame range for the majority of the 
sample. An example of the fit is shown in figure |A.1[ The 



10000 

A [A] 


Fig. A.l. Linear fit (red line) to the best-fit SED in the plane 
log{X),log{Fx). 

errors related to the Pmodeii i-e. the UV slope derived 
from the MAGPHYS SED, are computed from the lin¬ 
ear fit. 


Appendix A.2: fd from the observed photometry 

In order to test the validity of the estimate of Pmodei 
we computed the UV-slope also from the observed 
data, when is available the photometric coverage in 
the rest-frame range of interest. We interpolated the 
observed photometry in the rest-frame rang e A € [1200 — 
3500] A, following the method described in 


Nordon et al. 


(20131. The photometric UV-spectral slope is defined as: 


(A.l) 


_ -0.4(Mi - M 2 ) , 

log{Xi/X2) 

where Mi and M 2 are the AB magnitudes at wavelengths 

Ai=1600 A and A 2=2800 A. 

The rest frame 1600 A and 2800 A magnitudes (Migoo and 
M 2800 ) were estimated by interpolating between the avail¬ 
able photometric bands. To derive Migoo we used the filters 
between rest-frame A G [1200 — 2800] A and interpolated 
the flux at I 6 OOA (converted to AB magnitudes) by fitting 
a linear function between the observed photometric 
bands. To derive M 2800 we selected the filter s tha t ob¬ 
serve the rest-frame A S [1500 — 3500] A. Figure A .2 shows 
an example for the computation of Migoo and AI 2800 for a 
source that have the required photometric coverage. Fig¬ 
ure A.3| shows the case in which we haven’t the coverage 
in the rest-frame. The jdphot was computed for 13 sources. 
We computed the errors for Pphot by using a set of MC 
simulations: we randomly varied the observed photometry 
within the errorbars and we then computed the interpola¬ 
tion for each set of simulated values of the observed pho¬ 
tometry. The la uncertainty on jdphot was then estimated 
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from the width of the probability distribution function of 
the simulated values, assuming that this distribution has a 
Gaussian shape. The median error on f3phot was added 
in quadrature to the error of Pmodei thus the error 
for the UV-slope is cr|_. = + median 



value strongly dependent on the model to those constrained 
by the observed data. The correlation between Pmodei and 
Pphot is shown in figure A.4 the agreement betweeen the 
two estimates is quite good, considering also the poor statis¬ 
tic, and confirms the reliability of Pmodei derived by fitting 
the modeled SEDs. 



^ phot 


Fig. A.2. The plot displays the method for the computation of 
Pphot for the source 3213 (ID MUSIC). In the figure it is specified 
also the redshift of the source. The black open diamonds are 
the observed photometry, connected with a linear interpolation 
(black solid line). The green vertical lines highlight the rest- 
frame spectral range considered to compute Mi while the blue 
vertical lines mark the range for the computation of M 2 . The 
two dash-dot red lines mark respectively the positions of 1600 
A and 2800 A rest-frame. 


Fig. A.4. Comparison between Pmodei, derived from the SED 
fitting, and Pphot, computed by fitting the observed photometry, 
for the 15 sources with photometric coverage in the rest-frame 
range A £ [1200 — 3500]A. The blue line is the 1:1 line. The 
linear Pearson correlation coefficient is r=0.96 


A rest-frame fAl 

0 1000 2000 3000 4000 



Fig. A.3. The plot is analogous to the previous. In this case we 
do not have the photometric coverage to derive Pphot- 


Appendix A.3: jSphot vs jSmodei 

The UV spectral slope is model-dependent, since it is ob¬ 
tained from a fit to the SED: in order to verify the validity of 
our measurements we compared Pmodei to Pphot for the 13 
sources that have the required photometric coverage in the 
UV spectrum. The Spectral Energy Distributions were re¬ 
computed in this case excluding the photometric bands in 
the UV rest-frame spectral range, as we want to compare a 
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Appendix B: Main parameters of the sample 

Table [BT| summarizes the MUSIC ID, the coordinates, the 
redshift measured from the 3D-HST near-IR spectra, the 
observed Ha luminosity I-iHa,obs corrected for the aperture 
as explained in section [3?^ the infrared luminosity Ijir, the 
observed UV luminosity Ljjv and the stellar masses M* of 
the sample. 
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Fig. 11. Comparison between SFR/ij+r/v and SFRrfa when varying the dust correction for the Ha emission. The lower panel 
shows the case in which Lrra are corrected for dust attenuation by using the classical prescription of Calzetti {f=0.58) while in 
the upper panel the Ha luminosities are corrected with our dust correction {f=0.93). The black solid line is the 1:1 correlation 
line and the different colors mark the stellar masses, as indicated in the vertical color-bar. In the lowest part of each panel are also 
shown the median errors on the SFR measurements. 

Table 3. Different values of / obtained from literature. In the table are also specified the redshift ranges (third column of Tab. 

t , the types of sample and the methods implied to measure the differential extinction on nebular lines (fourth and fifth columns 
Tab.^ respectively). 


Author 

/ 

^ range 

Sample 

Method 


Calzetti et al. 

(2000 

1 0.44 (0.58) 

0.003-0.05 

starburst galaxies 

B aimer decrement 

1 

Kashino et al. 

|20i3 

1 0.83T0.10 

1.4-1.7 

sBzk galaxies 

Balmer decrement (stacked spectra) 



0.69T0.02 



Ha to UV SFRs 


Wuyts et al. 

20111 

0.44 

0-3 

Kg selected galaxies 

SFR indicators 


Price et al. ( 


0.55±0.16 

1.36-1.5 

MS star forming galaxies 

Balmer decrement (stacked spectra) 

1 

^annella et al. 

(2014 

1 0.58 

< 1 

UVJ selected galaxies 

comparison between Aj/y-M*/A//Q,-M* 

II 


0.77 

1 

1 

> 1 

II 

II 

\ 

tylentino et al. (201f 

)) 0.74±0.05 

2 

CL J1449+0856 

Balmer decrement (stacked spectra) 


hirb et al. {z 

'Dour 

1 

2 

rest-frame UV selected galaxies 

matching UV and Ha SFRs 


Reddy et al. 

20101 

1 

1.5-2.6 

Lyman Break Galaxies 

matching X-ray, UV and Ha SFRs 

1 Jiis wor 

k 

0.93T0.06 

0.7-1.5 

far-IR selected galaxies 

Ha to UV SFRs 
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Table B.l. List of the main parameter measured for the galaxies in the sample. The first column indicates the MUSIC-ID of the 
sources, the second and third columns are the coordinates, the 4*^ column reports the redshift derived from the 3D-HST near-IR 
spectrum, the 5*^ column is the observed Ha luminosity (i.e. not corrected for dust attenuation), the 6^^ column reports the 
infrared luminosity, the 7*^ column is the bolometric rest-frame UV luminosity (derived from Lieoo, as explained in Section 4.21 
and the last column is the stellar mass. For each of the quantities listed here it is also indicated the Icr uncert ainty . 
The quantities L/rj, huv and M* are derived from the SED fittin g w ith the MAGPHYS code (see Section ETm. 
The Ha luminosity is corrected for the contamination by [Nil] (cfr. Section [4.2|) and for the aperture (see Section [3. 2 [| 


MUSIC ID 

RA 

[decimal deg.] 

Dec 

[decimal deg.] 

Z3D-HST 

log(Li/Q_of,s)±l(T 

[erg/s] 

log(L7ij)Tlcr 

[Lol 

log(Li6oo)±le' 

[Lol 

log(M*)Tlcr 

[Mq] 

3213 

53.0714684 

-27.8809948 

0.99 

41.95±0.03 

11.46T0.07 

9.71T0.97 

10.06T0.10 

3636 

53.0941353 

-27.8755627 

0.79 

41.37T0.06 

10.97T0.06 

9.81T0.98 

10.21T0.07 

3698 

53.0734024 

-27.8745785 

1.10 

42.34±0.05 

11.82T0.03 

9.72T0.97 

11.30T0.06 

3827 

53.0760498 

-27.8736115 

0.74 

41.98±0.01 

11.20T0.01 

9.34T0.93 

10.32T0.01 

4491 

53.0893517 

-27.8644161 

0.83 

42.27±0.01 

10.89T0.09 

9.99T1.00 

9.90T0.25 

4738 

53.2219315 

-27.8594131 

0.74 

42.05T0.02 

10.82T0.05 

9.90T0.99 

10.46T0.09 

5142 

53.2202797 

-27.8540936 

1.00 

41.85T0.04 

11.72T0.02 

9.66T0.97 

11.10T0.07 

5302 

53.2209053 

-27.8513260 

0.74 

42.25T0.01 

11.14T0.04 

8.13T0.81 

10.76T0.09 

5334 

53.1965446 

-27.8516655 

1.22 

42.00±0.04 

11.74T0.01 

9.17T0.92 

11.16T0.08 

5429 

53.0707664 

-27.8506660 

0.99 

41.57±0.42 

11.36T0.02 

8.62T0.86 

9.84T0.06 

5812 

53.1155014 

-27.8446636 

1.09 

42.00±0.07 

11.82T0.02 

9.62T0.96 

11.08T0.06 

6272 

53.1077652 

-27.8388157 

1.11 

41.74T0.14 

11.40T0.01 

9.96T1.00 

10.33T0.01 

6344 

53.2149391 

-27.8382416 

0.85 

42.38T0.01 

11.24T0.01 

10.24T1.02 

10.02T0.08 

6722 

53.2181320 

-27.8336830 

0.85 

42.15±0.03 

11.24T0.01 

9.57T0.96 

10.31T0.10 

6983 

53.2155266 

-27.8308353 

0.74 

42.16±0.02 

11.33T0.01 

9.66T0.97 

10.23T0.08 

7224 

53.0898857 

-27.8282318 

1.10 

42.09±0.02 

11.17T0.04 

9.32T0.93 

10.11T0.17 

7335 

53.0936508 

-27.8263798 

0.88 

42.30T0.01 

11.93T0.01 

8.84T0.88 

11.32T0.09 

7400 

53.2023735 

-27.8262043 

1.11 

41.48T0.10 

12.03T0.01 

7.81T0.78 

11.05T0.01 

8003 

53.0892830 

-27.8171749 

0.75 

42.60T0.01 

10.82T0.04 

9.83T0.98 

9.84T0.06 

8090 

53.0667534 

-27.8166199 

1.41 

41.84±0.11 

12.11T0.01 

8.14T0.81 

11.31T0.02 

8231 

53.1498871 

-27.8140030 

1.31 

41.69±0.12 

11.09T0.04 

10.52T1.05 

10.34T0.20 

8737 

53.0681343 

-27.8066711 

1.30 

41.93±0.17 

11.41T0.01 

9.60T0.96 

lO.OlTO.Ol 

9269 

53.1198807 

-27.7987385 

1.38 

42.45±0.02 

11.59T0.01 

9.35T0.93 

10.85T0.08 

9302 

53.0341835 

-27.7977962 

0.84 

42.24T0.01 

11.45T0.02 

9.46T0.95 

10.58T0.12 

9310 

53.0292320 

-27.7981472 

0.84 

41.91T0.01 

10.92T0.05 

8.69T0.87 

10.45T0.09 

9452 

53.1567535 

-27.7956085 

1.11 

42.10T0.04 

11.22T0.03 

9.73T0.97 

10.31T0.07 

9687 

53.1841621 

-27.7926331 

0.73 

41.14T0.19 

11.16T0.02 

9.79T0.98 

9.80T0.05 

9779 

53.0268173 

-27.7913227 

1.02 

42.18T0.03 

11.71T0.02 

9.61T0.96 

11.18T0.09 

9903 

53.0906639 

-27.7901840 

0.99 

41.82T0.04 

11.15T0.05 

9.33T0.93 

10.29T0.13 

10015 

53.1095657 

-27.7882004 

0.99 

42.29T0.01 

10.92T0.09 

9.99T1.00 

9.89T0.05 

10263 

53.1765938 

-27.7854614 

1.31 

42.03T0.02 

11.45T0.01 

8.84T0.88 

11.15T0.04 

10295 

53.0440521 

-27.7850533 

1.04 

41.78T0.03 

11.40T0.03 

9.81T0.98 

10.17T0.16 

10617 

53.1352501 

-27.7816753 

1.43 

41.91T0.05 

10.54T0.10 

9.26T0.93 

10.32T0.06 

10669 

53.0406151 

-27.7803345 

0.96 

41.91T0.02 

10.98T0.08 

9.31T0.93 

10.04T0.10 

10674 

53.1790581 

-27.7805214 

1.03 

42.21T0.01 

11.73T0.04 

7.55T0.75 

11.04T0.10 

10681 

53.1288300 

-27.7804184 

1.17 

42.58T0.01 

11.29T0.02 

9.03T0.90 

10.67T0.08 

10730 

53.0497627 

-27.7790527 

1.22 

41.45T0.il 

11.44T0.02 

9.24T0.92 

10.91T0.05 

10859 

53.1452293 

-27.7779026 

1.09 

42.01T0.02 

11.56T0.02 

9.11T0.91 

10.44T0.07 

10972 

53.0584106 

-27.7761879 

0.83 

42.01T0.01 

11.44T0.02 

8.54T0.85 

10.26T0.06 

11121 

53.0370216 

-27.7747364 

1.03 

42.17T0.01 

10.07T0.07 

8.94T0.89 

10.37T0.10 

11346 

53.1279831 

-27.7714329 

1.31 

41.59T0.06 

11.21T0.04 

9.90T0.99 

10.71T0.06 

11381 

53.1059380 

-27.7714195 

0.89 

41.92T0.02 

11.46T0.02 

9.39T0.94 

10.60T0.10 

11658 

53.1169853 

-27.7683487 

1.11 

41.77T0.04 

11.37T0.03 

8.44T0.84 

10.69T0.07 

11708 

53.0656624 

-27.7678661 

1.53 

42.06T0.04 

11.59T0.04 

9.61T0.96 

10.30T0.10 

11794 

53.1513672 

-27.7666893 

0.89 

42.20T0.01 

11.08T0.01 

8.94T0.89 

10.23T0.01 

11814 

53.1812057 

-27.7656651 

1.22 

41.60T0.12 

11.37T0.06 

9.80T0.98 

10.44T0.06 

11958 

53.0285225 

-27.7639732 

0.84 

42.40T0.01 

10.95T0.04 

9.83T0.98 

9.96T0.07 

12297 

53.1636009 

-27.7589455 

1.09 

41.71T0.04 

11.29T0.02 

9.35T0.93 

11.11T0.06 

12349 

53.1119385 

-27.7578316 

0.83 

42.33T0.01 

11.15T0.02 

9.00T0.90 

10.87T0.07 
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Table B.2. List of the main parameter computed for the galaxies in the sample (continuation). 


MUSIC ID 

RA 

[deg] 

Dec 

[deg] 

Z3D-HST 

log(L//Q,^o&s)±l(T 

[erg/s] 

log(L/H)±lcr 

[L©] 

log(Lieoo)±l^ 

[L©] 

log(M*)±la 

[Mq] 

12401 

53.0729218 

-27.7578526 

1.04 

42.15±0.01 

11.04±0.01 

9.56±0.96 

10.30±0.01 

12468 

53.0228424 

-27.7570210 

1.19 

42.07±0.02 

11.63±0.02 

8.64±0.86 

10.95±0.07 

12660 

53.0804939 

-27.7538910 

0.73 

41.79±0.03 

10.78±0.07 

9.32±0.93 

9.75±0.06 

12667 

53.1952438 

-27.7537766 

0.84 

41.94±0.01 

11.26±0.08 

9.34±0.93 

10.56±0.05 

12781 

53.0962143 

-27.7525444 

1.01 

42.91±0.01 

10.39±0.06 

8.92±0.89 

10.36±0.07 

12815 

53.0361824 

-27.7522297 

1.29 

42.15±0.04 

11.92±0.01 

8.95±0.89 

11.22±0.04 

12967 

53.1073837 

-27.7498131 

0.83 

42.06±0.01 

10.99±0.02 

9.36±0.94 

9.64±0.06 

13192 

53.1616096 

-27.7469234 

0.73 

41.59±0.02 

11.72±0.01 

8.38±0.84 

10.91±0.01 

13194 

53.0687332 

-27.7469501 

0.97 

41.38±0.12 

11.44±0.02 

9.66±0.97 

10.54±0.16 

13382 

53.1463699 

-27.7443180 

1.03 

41.66±0.03 

11.29±0.02 

9.70±0.97 

10.18±0.03 

13505 

53.0615921 

-27.7425251 

0.73 

42.37±0.01 

11.03±0.13 

7.84±0.78 

10.37±0.05 

13540 

53.1073494 

-27.7419224 

0.89 

41.47±0.06 

11.48±0.02 

7.93±0.79 

11.03±0.18 

13552 

53.0940933 

-27.7405052 

0.73 

41.45±0.06 

11.49±0.01 

9.86±0.99 

11.11±0.03 

13617 

53.1501236 

-27.7399464 

1.05 

42.11±0.03 

11.30±0.04 

10.02±1.00 

10.12±0.09 

13664 

53.0837898 

-27.7395668 

1.22 

42.08±0.03 

11.51±0.02 

9.48±0.95 

10.93±0.09 

13694 

53.1752434 

-27.7392445 

1.15 

41.29±0.23 

11.10±0.07 

9.50±0.95 

10.31±0.09 

13844 

53.0494080 

-27.7370377 

1.28 

41.87±0.09 

11.43±0.06 

9.86±0.99 

10.51±0.06 

13915 

53.1443443 

-27.7355614 

1.26 

41.90±0.04 

11.11±0.02 

8.94±0.89 

10.80±0.08 

14519 

53.0166168 

-27.7270679 

0.96 

42.52±0.01 

11.57±0.02 

8.80±0.88 

11.23±0.08 

14585 

53.0393867 

-27.7266464 

1.02 

42.41±0.01 

10.91±0.01 

8.74±0.87 

10.65±0.01 

15143 

53.1693153 

-27.7189751 

1.22 

41.86±0.07 

11.39±0.01 

9.94±0.99 

10.65±0.01 

15279 

53.1508827 

-27.7161102 

0.96 

42.37±0.02 

11.37±0.02 

10.53±1.05 

10.13±0.08 

15568 

53.1889420 

-27.7137222 

1.19 

41.89±0.13 

11.49±0.07 

9.16±0.92 

10.83±0.09 

15711 

53.1235085 

-27.7118225 

0.65 

42.16±0.01 

11.38±0.02 

9.79±0.98 

10.27±0.11 

16046 

53.0786095 

-27.7074299 

1.08 

41.48±0.08 

11.13±0.05 

9.56±0.96 

10.09±0.07 

16930 

53.1391869 

-27.6941433 

1.04 

42.47±0.01 

11.45±0.02 

9.37±0.94 

10.96±0.05 

17069 

53.1626816 

-27.6923409 

0.73 

42.52±0.01 

10.39±0.09 

9.76±0.98 

9.55±0.05 

17450 

53.0881500 

-27.6851749 

1.12 

42.75±0.01 

11.04±0.04 

9.52±0.95 

11.01±0.07 

17593 

53.1090813 

-27.6835785 

1.04 

42.53±0.01 

11.04±0.08 

9.49±0.95 

9.18±0.09 

30257 

53.1665344 

-27.7037868 

0.82 

41.44±0.06 

10.93±0.07 

9.42±0.94 

9.69±0.04 
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